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We evaluate the factor k appearing in Langer's expression for the nucleation rate extended to the 
case of hadron bubbles forming in zero baryon number cooled quark-gluon plasma. We consider both 
the absence and presence of viscosity and show that viscous effects introduce only small changes in 
the value of k. 

PACS numbers: 24.85. +p, 12.38. Mh, 25.75.+r, 64.60. Qb 



I. INTRODUCTION 

There has been a great deal of recent activity related 
to the study of the quark-hadron phase transition. For 
the transition from quark matter to hadron matter there 
is interest in the rate of formation of critical size hadron 
bubbles, specifically bubbles of pion gas, in a supercooled 
vapor [1,2]. The nucleation rate is an essential ingredient 
of the complete treatment of a phase transition, as criti- 
cal size bubbles must first be formed in the quark vapor 
before the vapor can be converted into hadron material. 
Hadron bubbles initially form due to fluctuations of the 
energy density in the vapor. Those with radii, R, smaller 
than a critical size Ro collapse, while those of the critical 
size begin to grow exponentially. Langer and Turski [3] 
have shown that the nucleation rate, /, can be written as 
follows: 



2tt 



r2oexp(-AF/T). 



(1) 



Here k is related to the growth rate of the radius, R, 
of the bubble near the critical radius, Ro, by dR/dt « 
k{R — Ro). The factor is a statistical prefactor, and 
AF the difference in free energy of systems with, and 
without, a critical size drop present. This paper is con- 
cerned primarily with discussing the rate k. 

Csernai and Kapusta [1] have applied the Langer- 
Turski nucleation rate, (1), to the formation of hadron 
bubbles in a quark vapor. They consider both the hadron 
and quark materials as substances with zero baryon num- 
ber and treat them using relativistic formalism. One 
principal result of [1] is the suggestion that for baryon- 
free matter k is proportional to the viscosity coefficients 
of the quark vapor, and thus when these coefficients van- 
ish a hadron bubble does not form. 

We reexamine the relativistic evaluation of k in this 
paper, and show that viscosity is not necessary for the 
growth of hadron bubbles, and that its effects are of 
higher order than the energy flow terms in the growth 
process. The difference between our results and those 
of Ref. [1] seems to arise from the conceptual distinction 



between heat conduction and energy flow. The former is, 
indeed, not defined in baryon free systems i.e. systems 
with a chemical potential of zero as pointed out in [1] and 
other sources. Energy flow, however, is defined. A dif- 
ference between our results and Ref. [1] also arises from 
the technical difference in the treatment of the pressure 
gradients. 

The paper is organized as follows. In Section 2., we 
review the treatment of energy flow, viscosity and ther- 
mal conduction in the baryon, and baryon- free, cases 
of relativistic hydrodynamics. We consider the distinc- 
tion between heat conduction and energy flow and stress 
that the latter can occur without viscosity in the baryon- 
free case. In Section 3., we calculate the rate k in the 
absence of viscosity and baryon number by solving rela- 
tivistic energy and momentum transport equations. The 
energy and momentum relativistic equations are mathe- 
matically of the same form as the mass conservation and 
momentum nonrelativistic transport equations for a non- 
viscous vapor with perfect heat conduction. We compare 
the relativistic result with the solution for k in such a 
nonrelativistic case as obtained in [3] , and show that the 
forms are similar. The method of solution of the trans- 
port equations used in this paper differs from that used 
in Ref. [3]. We compare the two methods, both of which 
can be used in either the relativistic and nonrelativistic 
cases. In Section 4, we reexamine the calculation of the 
rate k in [1] , focusing on the relationship between the free 
energy and the gradient of pressure, which is used in the 
momentum transport equation. In particular, we suggest 
that a source for the difference between the conclusion of 
this work and those of Ref. [1] , is linked to a difference in 
treatments of the pressure gradients. Finally, in Section 
5, we examine the relativistic transport equations in the 
presence of viscosity and suggest that viscous effects are 
of high order in the case of small viscosity coefficients. 
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II. ENERGY VERSUS HEAT FLOW IN BARYON 
AND BARYON-FREE RELATIVISTIC 
HYDRODYNAMICS 

We first review the transport relationships for a system 
with nonzero baryon number. This is a system for which 
the baryon chemical potential is also nonzero. Under 
this condition there is a net number of either particles 
or antiparticles present . Thus one can define a number 
density, n, and obtain the conservation of net particle (or 
antiparticle) number in the absence of dissipative effects: 



and for i ^ Q 



d{nu^) 



= 0. 



(2) 



Here represents the 4- velocity (divided by the velocity 
of light) of the net particles. The energy density and 
flow are given by the stress tensor (also in the absence of 
dissipative effects): 



(3) 



which provides energy and momentum conservation 
through the relationship: 



dx^ 



0. 



(4) 



Here g^^ is the metric, e, the energy density, p, pressure, 
and Uo = \l \f{\ — ?;^/c^) = 7. The combination e + p is 
the enthalpy density, w. Both the energy flow and parti- 
cle flow are related to m*. A point we wish to emphasize 
is that the term wu^u^ represents the energy flux in the 
absence of viscosity. The enthalpy density depends on 
the material (particles and antiparticles) in the respec- 
tive volume. 

For the case of nonzero dissipation, an additional term 
T*" containing viscosity coefficients (see [5]) is added to 



in (|). The new tensor still satisfies (|j). However 
the particle conservation relationship (0) is replaced by: 



= 



(5) 



where is proportional to the coefficient of thermal con- 
duction [5] . In the relativistic dissipative case it is conve- 
nient to associate with the velocity of energy flow, and 
to note that an energy flux necessarily involves a mass 
flux [5]. 

Next we examine the case of zero baryon number, i.e. 
the case of zero chemical potential and n = 0. In this case 
there is no net particle number and the thermal conduc- 
tivity is not deflned, because there is no net baryon den- 
sity with respect to which energy can be conducted. The 
energy flow, however, need not be zero. The energy flow 
is explicitly contained in (^. In the low velocity limit, 
i.e. to first order in v and in the absence of viscosity, 
equation (Eh gives: 



dtwv 



-Vp. 



(7) 



Here, the velocity of light has been set equal to 1. We 
wish to emphasize that the energy flow, wv, given by (|^), 
does not vanish, even though there is no heat conduction. 

In the presence of viscosity, terms of second order in u 
appear in the energy equation (^) , while a term linear in 
u appears in (^, the momentum equation. This means 
that the viscosity terms are relatively unimportant in the 
energy transport when u is small. The momentum equa- 
tion, however, indicates that viscosity influences the time 
evolution of u. Thus viscosity can serve to disrupt the 
energy flow and generate entropy but cannot be a driving 
mechanism for energy removal. 

The authors of Ref. [1] also list equation (^) with the 
high order viscosity terms dropped. This explicitly shows 
that energy can change in time in the absence of viscos- 
ity. In Ref. [1], however, a different equation with the 
dominant term oc wu^ missing, and with the second order 
viscosity term present, is also given and used to describe 
the rate of energy change. 



. dR 4 duR 



(8) 



dte — —V • [wiP) 



(6) 



The terms 77 and C are the shear and bulk viscosity co- 
efficients respectively, and the left hand side represents 
the energy fiux density (^energy per unit area per time). 

Equations (^), and ^ are inconsistent. Equation 
(^)suggests that the energy flow is provided solely by 
viscous effects. It further suggest that, in the absence of 
such effects, the hadron bubble could not grow because 
the flow of the latent heat generated by the phase transi- 
tion is forbidden. The fact that a v is present, however, 
means there is energy flow oc ww*, regardless of viscous 
effects, and the fact that heat conduction can not be 
deflned. Physically one may suggest that the energy is 
being carried away by quarks and antiquarks i.e. by con- 
vection to list a possible mechanism. 



III. SOLUTION OF RELATIVISTIC FLUID 
EQUATIONS FOR NONVISCOUS CASE 

We now consider solving explicitly for the rate k in the 
absence of both viscosity and baryon number. This exer- 
cise, explicitly shows that k does not vanish. In Section 
5, we further suggest that the addition of viscous effects 
leads to small corrections in k. 

In the nonviscous case and in the limit of small veloc- 
ities the relativistic energy and momentum differential 
equations become [5] (^ and (^. These equations have 
the same form whether the baryon number is zero or not. 

The method of solution we use here was suggested by 
C. Goebel. The differential equations are solved in the 
quark vapor region only and are then related to the inte- 
rior hadron material through boundary conditions. One 
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begins with the hnearized relativistic equations in the 
quark region for a velocity field with spherical symme- 
try: 
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dte + —Wydrir v) ^ 
WydtV + drP — 0. 



(9) 



The boundary conditions or Kotchine conditions, across 
the interface are found by equating T''' in the two regions 
in a frame in which the surface of discontinuity is at rest 
and also T" in the two regions with a surface tension 
term added. The results are: 



Ph 

Wy{Vy - 



2(7 

R) = Whivh - R)- 



'Pv 



(10) 



Use has been made of the Galilean transformation as 
V << I and only terms linear in v retained. The v's 
above represent velocities measured relative to the bub- 
ble center. Here the subscript v refers to the quark vapor 
and h to the hadron vapor. R is the velocity of the bub- 
ble surface. For simplicity the velocity inside the hadron 
bubble, Vy is taken to be zero for the following solution. 
This condition is not an essential one (see Appendix A). 
A third boundary condition is needed to give information 
about the combustion rate. We postulate it as: 



wyivy-R)^cA{sn-sn). 



(11) 



Here C is a constant related to the efficiency of energy 
flow across the quark-hadron interface. For illustration 
purposes we take: 



STh — STv 



(12) 



which is equivalent to C — > cxd. In a future work we 
examine the effect of finite C. Since both e and p are 
functions of the temperature, T, alone, the variations in 
e and p, i.e., e and p, satisfy 



P 



where 



(13) 



(14) 



the prime denoting the T derivative. It can be shown 
that c is the velocity of sound in the vapor. Next one can 
considers a variation of pressure away from the stationary 
solution, p. Then: 



Pv -Ph 



2aR 



(15) 



Here Ro is the critical size and R the variation from this 
size. The quantities are all evaluated at Ro in (|lO|). 

Given the following forms for vapor and hadron pres- 
sures: 



Py = ayT'^ - B 
377r2 

Ph = a/iT^ 



cth = 



3^ 
90 ' 



finds: 



Pv 
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at 



Rl 



(16) 



(17) 



if one takes 5Ty = STh. This particular choice provides 
for changes in pressure in both the hadron bubble and 
quark vapor regions as the bubble grows. Differentiating 
( p3| ) with respect to time and using (17) one obtains a 
relation which holds at Ro : 
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where 



and 
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T^oi^uj) 



WyV 



(18) 



(19) 



(20) 



if Vh is taken to zero as a first approximation (see Ap- 
pendix A) . In the vapor region (Hq) becomes a wave equa- 
tion upon differentiation with respect to time and sub- 
stitution of the second relation of (^ differentiated with 
respect to r : 



Explicitly using: 



leads to: 



/'(r - ct) 



r{rf'-f)- 



(21) 



(22) 



(23) 



At this stage / is an arbitrary function and /' is its 
derivative with respect to r — ct. Using ( |l8|) one obtains: 



/" + y/'-^/ = o. 

Ho 

Solving with the form f(x) — exp(/ij:) one obtains: 



Ro 



)• 



(24) 



(25) 



The upper sign is associated with the growing mode. 
Finally, the rate k can be associated with — c/i, so that 
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In the limit of small y one obtains: 



(26) 



(27) 



where N is constant one can obtain p from 



For the purpose of comparison, we next review the so- 
lution provided in Ref. [3] for a nonrelativistic system 
with nonzero chemical potential. That work provides a 
K for a drop (e.g. water) growing in a nonviscous vapor 
with perfect heat conduction. ( Note: There is an error 
in the second part of [3] which involves thermal conduc- 
tion. The error is corrected in [7] and is unrelated to the 
first evaluation of n in [3] which ignores changes in tem- 
perature and hence also the specific thermal conduction 
terms.) The relativistic equations (^) solved above are 
of the same form as the nonrelativistic mass conserva- 
tion and momentum equations used in the first part of 
[3]. Thus one can compare the methods of solution and 
the results for the nonrelativistic treatment of [3] with 
the relativistic case. 

The nonrelativistic equations of [3], expanded about a 
stationary solution, are: 



dtv — —V • (nu) 
ndtu — — Vp. 



(28) 



Here n is the stationary density and v and u the changes 
in density and velocity from the stationary results. It is 
further assumed that the time dependence of all quanti- 
ties is ex-p{Kt). The term Vp is obtained in [3] from the 
following expression for the free energy: 



Fk - 
Fi - 

together with: 



F = Fk + Fi 
1 



d\[^K{Vnr + f{n)] 



1 ^SFi 
-p = V — 

n on 



(29) 



(30) 



where n is evaluated at n + iy. Equation ( |30| ) is consistent 
with the thcrmodynamical relation: 



where 



p^fiN-f 



dF 



^~ dN 

holding T and V constant. Using 

N 

F = fiy)V 



(31) 



(32) 



(33) 



P = 



dF 
"dV 



holding N and T constant. Thus 
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One may then obtain: 



1 



dr 



dtu = — v(-/-i:v^ 



dn 



(34) 



(35) 



(36) 



It is interesting to note that although the expressions 
for the pressure in terms of / are different for the zero- 
baryon number relativistic and nonzero baryon number 
nonrelativistic cases the change in pressure similar. For 
the two cases one has: 



P 
P 



c?e 



(37) 



respectively. In both cases cf denotes the speed of sound 
in the vapor. For the nonrelativisitic case: 



dp 
dn 

dn^ 



(38) 



The first equation of ( psD is of the same form as the 
first equation of (|^) but the two have different origins. 
The latter is the result of energy conservation which fol- 
lows from (^) involving the energy-stress tensor. The 
former is a result of particle conservation. There exists 
in relativistic fluid mechanics for systems with nonzero 
chemical potentials an equation similar to the first of 
( p8| ) which describes the conservation of a net number of 
either particles or antiparticles. Such an equation only 
holds for nonzero baryon number and is not used in this 
paper which deals with the zero baryon number case only. 
As a result it is better not to look for a one-to-one corre- 
spondence between the variables of the relativisitic and 
nonrelativistic equations even though ( p8| ) and (||) arc of 
similar form with n replacing w. 

One can proceed to solve the wave equation resulting 
from a manipulation of ( ^8|) in all space instead of solving 
just in the exterior region and using the Kotchine bound- 
ary conditions as done in the method of section 3. The 
wave equation is straightforward to solve in the exterior 
and interior regions of the bubble. In the interface region, 
however, the solution is more difficult as d"^ f /dn? varies 
and Langer and Turski [3] solve it approximately using 
the assumption of small n. Finally one may integrate v 
over all space and set this value equal to zero because 
the total amount of material is conserved. The integral 
in the interior contributes little and is ignored in [3] . The 
result of [3] then gives: 
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mRHAny 



(39) 



Here riy is the vapor density. It is interesting to note 
that K vanishes with vanishing Uy. The factor ny/An 
enters the problem in the following way. In the interface 
region the solution of rv taken from [3] is equal to: 



{ — ) 



R^An dn 
2a 'dr 



at r near R. Thus: 
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R'^An dr 



.dn 



(40) 



(41) 



The LHS is Spy which is equivalent to py and the factor 
Uyj An is related to the expression for the pressure of the 
vapor in equilibrium with a drop in this case namely: 



Py {total) = Po + Spy 

2any 

= Po 



RoAn 



(42) 



where Po is the pressure of a system with a large amount 
of liquid in equilibrium with vapor. The term 5py de- 
scribes the change from this pressure due to the effects 
of the surface of the drop. The procedure of solving in 
the interface region simply ensures that the outside and 
inside pressures are related by: 



2a 

"r ■ 



(43) 



One may now attempt to apply the above ideas to the 
relativistic equations (pP). These lead to: 



-V^p. 



(44) 



In the nonrclativistic and zcro-baryon number relativistic 
problem Vp is equal to a constant multiplied by v and e 
respectively. Furthermore the same Kotchine equations 
( p^ ) with w's replaced with n's apply to the nonrclativis- 
tic problem. As a consequence the low k relativistic zero- 
baryon result must be the same as the nonrclativistic one 
with n's replaced with w's or namely: 



, 2aWy 
-RHAwj 



(45) 



This is identical to 



IV. FREE ENERGY AND PRESSURE 
GRADIENTS 



of with an additional term linear in viscosity coeffi- 
cients to describe viscous effects. The result of [1] that k 
is proportional to the viscosity coefficients follows from 
the replacement of Vp in (^) with a particular expres- 
sion involving the free energy of the bubble plus vapor 
system. We examine the relationship between pressure 
and free energy more closely to see how it leads to the 
conclusions of [1] concerning the role of viscosity. 
The free energy is taken to have the form B : 



F = Fk+Fi 



K 



— I d^rwv^ 



Fi^ J (Ve) 



/(e)]. 



(46) 



Here X is a constant related to a, the surface tension. 
The free energy is of the same form as that used in non- 
rclativistic work [3] but is taken as a function of e, the 
energy density, as opposed to n, the number density. In 
the baryon-free case there is no number density. In [1] 
(equations 17-25) the pressure gradient term is associated 
with: 



de[r 



= -K{V^e)Ve + 

oe 



(47) 



It is important to note that the 'Vp' is not simply a pres- 
sure but a combination of a pressure and a force term. 
Both of these terms are found in the Euler or Navier- 
Stokes equations. Such an identification must be made 
for the following reason. The pressure on the inside dif- 
fers from that on the outside so there is necessarily a 
pressure gradient in the interface. If there were no force 
term then the Euler equation would require a changing 
fluid velocity in the equilibrium situation which is un- 
physical. The force term is associated with the KiVe)^ 
term which is related to the surface tension by: 



47ri?V = K 



d^7 



(48) 



This is the term needed to balance the differing pressures. 
Identical considerations hold in the nonrclativistic case, 
with n taking the place of e. 

In a region of e(r) varying slowly equation(|47|) suggests 
the identification: 



9/ 
de 



Ve = V/ = -Vp 



(49) 



and 



The work of Ref. [1] suggests that energy is removed 
during bubble growth by viscous effects alone (H). This 
result is already suggested by the authors' treatment in 
[1] of the relativistic momentum equation, i.e. the second 



P 



(50) 



The term proportional to K has been dropped in this 
region. We also assume this form. 
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Let us examine the approach used in [1] . There one lin- 
earizes the relativistic momentum equation for smah vari- 
ations, e for the energy and v for the vefocity, about the 
stationary solutions of e and respectively for the hadron 
bubble in equilibrium with the quark vapor. Specifically: 



oe ae^ 



The assumption: 



oe 



(51) 



(52) 



is used. This is a consequence of the minimisation of ( [iq ) 
with respect to e to obtain a profile for the stationary 
solution e. The stationary solution of [1] is described by 
the solution of: 



dFi 

de 



= 



5/ 

de 



(53) 



where a fourth order polynomial in e is chosen for /. This 
approach follows nonrelativistic work where n is used in 
place of e. In the nonrelativistic case, however, one con- 
strains the number of particles present and so does not 



The energy and momentum fluid equations in the pres- 
ence of viscosity, keeping only first order terms in v the 
velocity, are: 



dte = -V • WV + 0{v^) 
4 

wdtv = -Vp + V[(-77 + C)V • v\ 



(58) 



These are the same as those of [1]. For small viscosity 
coefficients and small k, it appears that the viscosity term 
is small compared to other terms as it involves two gradi- 
ents. Here we assume a non- vanishing pressure gradient. 
From the equations with viscosity absent one can see that 
gradients of v are smaller than e by a factor of k, and 
K is assumed to be small. Next we consider the effect of 
viscosity on the boundary conditions. In the absence of 
viscosity the boundary conditions ( |l0[) can be expressed, 
in a frame moving with the interface, by 



In the presence of viscosity: 



pg- 



(59) 



(60) 



dl 

dn 



= 0. 



With the use of (|5|) equation (64) of [1] 



(54) 



dt{wv) = (Ve)[-K{V^) + /'> + V[(C + -??)V • v\ (55) 

is obtained. One next makes the association dtv — kv. 
Evaluating (^5|) in the quark vapor region where e is con- 
stant causes the first term in ( |55| ) to vanish and suggests 
K (X viscosity coefficients. 

Following this approach the velocity would not change 
in time and the bubble would not grow in the absence of 
viscosity. We evaluate (M) differently. The term: 



dl 

de 



(56) 



is not taken to be zero but rather c^. This does not 
cause Vp to vanish and so bubble growth occurs with or 
without viscosity. 



V. EQUATIONS WITH VISCOSITY PRESENT 

In the presence of viscosity a term t**^ is added to the 
tensor T'*^ (^) where: 



dx^ ~^ dx^ ~^ dx'- ~^ dx^ ^ 



to first order in velocity where only first order velocity 
terms in r are retained. 

In a spherical problem one has 



[p+i^v + Oy -iv-R) 

4?7(w - R) 



(61) 



Thus p in the first of (^0|) is replaced with the RHS of 
( |6l| ) while the second boundary condition of ( |lO| ) remains 
unchanged to first order in velocity. Here v is the velocity 
measured in a frame in which the bubble center is at rest. 
A Galilean transformation has been performed to relate 
this frame with that moving with the interface because 
all velocities are small compared to that of light. 

We first consider the effect of the viscosity-modified 
transport equations. One can solve ( pq) explicitly fol- 
lowing the method of section 3, using (Oh) and replac- 
ing — V • w in ( ^8| ) with dte/w. The equation in the vapor 
becomes: 



d^e = c V e-l-V — dte. 



Por constant in space one obtains: 



e = — exp(— -I- nt) 
r 



(62) 



(63) 



(57) 

where A is a constant and k and q are related by: 
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2 2 

c q 



2 

— Kq . 



(64) 



Thus in the vapor region one obtains the same result as 
( ^ ) but with c^, the speed of sound, replaced with: 



ih + ^ 



)■ 



(65) 



Due to the presence of a k factor the second term on the 
RHS of (^4|)is small compared to the first. This suggests 
that viscosity is a small effect. Solving the first of (||) 
which is the same for both the nonviscous and viscous 
situations an expression for v: 



H .,q 1 X , ^ 

V = — + ^)cxp(-gr + nt) 
wq'^ r r'^ 



(66) 



is found. 

To see the explicit effect of viscosity on k a relation- 
ship incorporating the boundary conditions is needed. 
Eq. (|6^ ) simply provides a relationship between two un- 
knowns q and k; an additional equation is needed relating 
the two. To obtain such an expression ( p|) a nd (|6l| ) are 
combined in the same manner as (^ and (|10|) in secton 3 
to obtain an expression analagous to (|l^). The addition 
of viscosity modifies (p7|). To see how this modification 
arises we consider the derivation of (^^ in section 3. To 
obtain the result ( p7| ) one solves for 6Th = ST^ = ST in 
the first equation of {nw by setting: 



(67) 



Py (X aySTy 



and taking a similar expression for ph- One then sub- 
stitutes the value for ST in ( |6^ ) to obtain py. In the 
presence of viscosity a similar method is employed but 
now UhjUy^R in ( |6l| ) must be replaced with expressions 
involving ST. This allows one to solve for ST and thus 
Py . Terms with V • v can be replaced with terms involving 
KB using the first of (|58|). Terms involving V • R are pro- 
portional to KCy because the second of ( |To| ) relates R to 
Vy . Terms involving Vh are dropped to be consistent with 
section 3 (see Appendix A) . Finally one must consider 
terms with R alone. Given ( |63| ) and ( |66[ ) it is possible to 
relate Vy and Cy at the interface to obtain: 



K . 1 , _ 

— (9+ 

q'^Wy R 



(68) 



The use of (^ij) essentially leads to the replacement of 
the factor: 



Wy 

Aw 



(69) 



appearing in ( pT[ ) with the factor 

.r , SfyK ifqH ifqH 

X = Wy/[Wh - Wyil + — 1 

Aw Aw Wq 

Uc^N^ , K I ... 

+ [VhWy - rivWh)(— + — ))J. 

RAwWy K cN Ro 



(70) 



Here 



fh — 2^7/1 + Oi 



N 



(1 



ih + o 



(71) 



Next Eq. ( [T^ ) is differentiated with respect to time to 
obtain an expression for py . The term p (|l3|) is replaced 

with an expression involving R which is in turn replaced 
with an expression in the velocity v. The term e is written 
in terms of the gradient of the velocity. One then obtains: 



—Wydr-r V — YwyV 



where: 



Y = yx- 



Aw 



(72) 



(73) 



y being given in (p^). The velocity in the hadron material 
is again neglected. 

Substituting into (|7^) yields: 



(74) 



Ho 



which may be compared with (^6|), the result in the ab- 
sence of viscosity. Eq. (^4|) must be solved with q ex- 
pressed in terms of k namely: 



(75) 



The modification due to viscosity resulting from the 
change in the transport equations (|5^) is present in q 
while that due to the change in the boundary conditions 
(|6l| ) is contained in the replacement of y with Y. In 
general (74) must be solved numerically for k as it is of 




higher than third order in k. An approximate solution is 
obtained if the q term can be neglected compared to 



(76) 



with Y and N being evaluated with k replaced with Kq, 
the result in the absence of viscosity (p5|). The solution 
for K in the absence of viscosity is of the same form as 
([TG]) except N is replaced with 1 and Y with y. 



VI. CONCLUSION 

We have examined the growth of a hadron bubble in 
quark vapor and suggest that the physical process is dif- 
ferent from that presented in [1]. Specifically we suggest 
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that the formahsm of relativistic hydrodynamics permits 
energy flow in the absence of heat conduction and vis- 
cosity. We have solved the relativistic fluid equations for 
a system with zero baryon number (no thermal conduc- 
tion) and no viscosity and have obtained a nonzero value 
for K. We have compared the relativistic fluid equations 
and resulting k with the nonrelativistic equations and k 
presented in [3]. The method used in [3] is approximate 
in that it only holds in the limit of small k and differs 
from that presented in this paper which provides an ex- 
act result for k. Finally we have solved the relativistic 
fluid equations for zero baryon number in the presence 
of viscosity. We suggest that viscous effects introduce 
only a small correction in the nonviscous value of k and 
provide an expression for the correction. 

In this paper we have used a specific set of boundary 
conditions linking the quark and hadron regions. Investi- 
gation of different sets of boundary conditions and their 
effect on k is in progress. In addition we consider the cal- 
culation of K for a medium with nonzero baryon number 
and thus heat conduction. 
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K = qch 



(A3) 



where Ch is the speed of sound in the hadron material 
and equals that of the quark medium. The solution in 
the quark region is of the form: 



Cg — Cq 



cxp[Kt — q(r — Ro)] 



(A4) 



Thus one may obtain the relationship: 



eh 



CLgSTg 



(A5) 



at the interface surface. Finally one may solve for k by 
evaluation of: 



R^Aw 



{WqVq - WhVh) 



(A6) 



with (A3). All quantities are evaluated at the surface of 
the bubble. The resulting expression for k is complicated 
due to the presence of the terms coshqR and sinhgi?. 
One may, however, consider the limit of small q which is 
equivalent to small k. In that case: 



Vh oc ~kR. 



(A7) 



APPENDIX A: A 



A solution for k has been obtained under the assump- 
tion that 6Th — STq at the surface of the bubble and 
that the velocity in the hadron material is zero. One 
may solve for k without the latter restriction to find a 
more general result which in the limit of small k reduces 
to the zero velocity solution. For small k the velocity on 
the inside is much less than that on the outside and so 
can be ignored. One may note that the authors of [3], 
who consider the problem of nonrelativistic drop forming 
in a vapor, also disregard the velocity inside the drop in 
the limit of small k. 

To solve the general problem with velocity retained in 
both the bubble and vapor regions the hydrodynamical 
equations must be solved in both regions and use made 
of the full boundary condition: 



WqVq - WqR = Wh^h ~ WhR- 

The solutions in the inner region are: 



(Al) 



Here the cosh and sinh terms have been expanded to sec- 
ond and third powers of q respectively. It is possible to 
see that vt is much smaller than the expression for Vq 
in the case of small q and so Vh = as an approxima- 
tion. One can procee d to solve for n with the next order 
approximation using (A7) instead of Vh = 0. Then one 
must solve the quadratic: 



2 2 

c q 



2aWa 



i?2(Azi;)2 



(9 



1 
R 



CLqi 



) 



to obtain: 



r 2 , / 4 2 4yc4 ahR, 



(A8) 



(A9) 



Here y is given by (|l9|). In the hmit y « ^ one obtains 
the results of section 3. 



APPENDIX B: B 



eh = Ch 



sinh(qr) e*^* 
sinh(qRo) r 

Ch K 



, , , sinh(qr) , e''' 

= • w p ^ 2 (-<ZCOsh(qr) + ^) — . (A2) 

smh(qKo) Whq r r 

In addition there is the relation: 



A method to solve for k using differential equations 
in a metastable and stable phase region together with 
Kotchine conditions relating the two phase at the bub- 
ble surface is presented in section 3. This method holds 
a, the surface tension, fixed as the bubble grows. An 
alternative method, which yields identical results, is con- 
sidered in this appendix. It is based on a modification of 
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Langer and Turski's approach to the nonrelativistic prob- 
lem. In the Langer and Turski approach, as outhned in 
section 3, one solves a wave equation for n in the vapor 
and interior drop regions and obtains solutions for n and 
V in these two regions. These solutions are functions of 
K. Next an approximate solution for n is obtained in the 
interface region assuming small values of k. Finally the 
relation: 



drn{r) = 0, 



(Bl) 



namely particle conservation is used to obtain a value for 
K. In [3] only low order terms in k are retained after 
integration. If one, however, evaluates the expression 
exactly and considers the manipulations in the interface 
to also be exact one obtains the same results as in section 
3. Specifically, in the nonrelativistic problem one has: 



dtfi 



-V • nv. 



The condition (Bl) is equivalent to: 



drdth — 

for the case of interest namely: 
dtfi = nh. 



(B2) 



(B3) 



(B4) 



The volume of all space can now be separated into three 
regions. Two of the regions are the inner drop and outer 
vapor regions and the third is the vicinity of the surface 
in which the density is changing. Using (B3) and Gauss' 
law one may make the identification: 



dtfi = inter facepiece + n„z;„ - UoutVout- 



The surface piece is equivalent to: 



{An)R. 



(B5) 



(B6) 



This describes the change in the interface as equivalent 
to a drop with a sharp surface growing by the amount R. 
The inner and outer areas defining the interface are taken 
as equal. An expression for R can be obtained from: 



Pvapor 



2any 



and 



Pv 



C?fl. 



(B7) 



(B8) 



Here the Kotchine condition relating pressures on the 
inside and outside regions has been used. Langer and 
Turski in [3] obtain a similar result by approximately 
solving the transport equation: 



q2 £ 



(B9) 



in the region in which the density is changing. The result 
of [ 3] k eeps only terms to first order in n but is equivalent 
to ( p^ and (B7). Thus one finally obtains: 



— riinVin l^out'^out 4^ 



c^hR^An 
2an„ 



(BIO) 



which is identical to the result of section 3 obtained using 
a method which utilised the Kotchine expressions and n 
and V solutions together with: 



P 



(Bll) 



Identical arguments hold for the relativistic problem. For 
the relativistic problem to first order in the energy flow 
velocity over the speed of light one has: 



dte 



-V • wv. 



(B12) 



This equation holds with or without viscosity for both 
baryon and baryon-free materials. Separating the volume 
again into three regions and using Gauss' law leads to the 
energy conservation condition: 



WhVh - WqVq 



= AwR. 



(B13) 
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